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ABSTRACT 

SDSS J1356+1026 is a pair of interacting galaxies at redshift z = 0.123 that hosts a luminous obscured 
quasar in its northern nucleus. Here we present two long-slit Magellan LDSS-3 spectra that reveal a pair of 
symmetric ~ 10 kpc-size outflows emerging from this nucleus, with observed expansion velocities of ^250 
km s" 1 in projection. We present a kinematic model of these outflows and argue that the deprojected physical 
velocities of expansion are likely ~ 1000 km s" 1 and that the kinetic energy of the expanding shells is likely 
1044-45 er g s -i^ .yyjj.jj an a b so i u te minimum of > 10 42 erg s" 1 . Although a radio counterpart is detected at 
1 .4GHz, it is faint enough that the quasar is considered to be radio-quiet by all standard criteria, and there is no 
evidence of extended emission due to radio lobes, whether aged or continuously powered by an on-going jet. 
We argue that the on-going star formation is probably insufficient to power the observed energetic outflow and 
that SDSS J1356+1026 is a good case for radio-quiet quasar feedback. In further support of this hypothesis, 
polarimetric observations show that the direction of quasar illumination is coincident with the direction of the 
outflow. 



1. INTRODUCTION 

The discovery of a tight relationship between the masses of 
black holes (BHs) in nearby galaxies and the velocity disper- 
sions and masses of their ste l lar populations ( Magorria n et all 
[l998t lGebhardtetaHl2000t iFerrarese & Merrittl 1 120001) sug- 
gests that the active phase of the supermassive BH evolution 
has profound effects on the formation of the galaxy. Indeed, 
there is direct evidence for such 'feedback' by radio-loud ac- 
tive galactic nuclei (AGNs) onto the large-scale gaseous en- 
viron ment of their host galaxies, both at low and a t high red- 
shifts dNesvadba et alJEOOH lFu"& Stocktonll2O09h . In these 
objects, ionized gas is outflowing with velocities above the 
escape speed from the galaxy, often along the direction of 
the radio jet. This also appea rs to be the mode in which 
AGNs heat up galaxy clusters (McNa mara & Nulsenl 120071) . 
However, as only a minority (^10%) of active BHs are radio- 
loud, this mechanism may not be sufficient to explain the col- 
ors and the luminos ity function of massive red galaxies (e.g., 
Sprin geTet al.l l2005). and it is as yet unclear if the remaining 
^90% of AGNs are affecting the gas, and thus the star forma- 
tion, in the host galaxies in a similar fashion. 

There is strong evidence that star formation driven winds 
play an important role in galaxy evolution. Their imprint is 
apparent in the mass- metallicity relation at low redshift (e.g., 
iTremonti et al.ll2004t) . and they are required i n galaxy sim- 
ulations to reproduce the observations (e.g., ISpringel et al.1 
12005b iDave et all 120081) . Most importantly, winds driven by 
star f ormation are clearly seen in low-redshift dwarf galaxies 
(e.g., [Martin 19 98b, in local ultra- l uminous infrared galaxies 
(ULIRGs; e.g. iRupke et al.ll2005al) . in galaxies at z ~ 1 (e.g., 
IWeiner et al. 2009) and high-redshift star-forming galaxies as 



well (e.g., iShaplev et all I2003T) . See iVeilleux et al.1 d2005l) 
for a recent review of low redshift observations. What has 
been much harder to pin down is th e role of radio-quiet 
BH accretion in drivi n g winds (e g., I Tremonti et al.] 1 2007 
Alexander etal] 120101: iKrugetalJEOlOb ICoil etaTT 12011 



Hainline et al. 201 1; Rupke & Veilleux 201 1). 



One class of BHs that may be showing radio-quiet feed- 
back processes in action are broad absorption line quasars, in 
which high-velocity outflows ap pear to be driven by the ra- 
diation of the quasar itself (e.g.. lProga et al.l 2000). Because 
the distance of the outflowing material from the BH and its 
covering factor are difficult to determine, the kinetic ener- 
gies of these outflows are unknown in most cases. In some 
sourc es, it may be a substantial fraction of the accretion en- 
ergy (ICrenshaw et al.l2003blMoe et al.l2009t) . and the detailed 
photo-ionization models place these outflows at a consider- 
able dis tance (several kpc) from the center of the galaxy (but 
see also Faucher-Gi guere et al.ll201 lb . These observations are 
promising, but so far have been possible only for a handful 
of sources. We seek unambiguous evidence for galaxy scale 
winds. 

A handful of strong cases have been made for radia- 
tively driven AGN outflows on gal a xy-wide scales at high 
redshift (e.g. I Tremonti eTall 12007b I Alexander etal] 120101: 
lHainline et al.l 1201 lb . Locally, new exciting observations 
of molecular lines, both with ground-based in terferometers 
dAlatalo et al.1 1201 ll) and with Herschel (e.g., iFischer et al.1 
1201 Oh are revealing large molecula r outflows t hat are plau- 
sibly AGN-driven. In particular, ISturm etal] (1201 ll) make 
a strong case that the molecular outflow velocity in local 
ultra-luminous infrared galaxies (ULIRGs) is strongly cor- 
related with the strength of the active nucleus, such that 
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more luminous A GNs drive faster outflows (although see also 
iKrug et alJl2OT0h . 

Perhaps the best low-redshift case study of large-scale out- 
flows is of Mrk 231. Mrk 23 1 is the most luminous lo- 
cal ULIRG (Sa nders et alJI 1 988h . the nearest broad emission- 
line AGN, a rare low-ionization broad absorption-line quasar, 
and contains a radio jet. This galaxy also hosts a massive 
young starburst that likely dominates th e bolometric output 
of the galaxy (e.g., lLon sdale et al. 2003). The galaxy has 
been studied in detail at e very conceivable wavelength (see, 
e.g., iFeruglio et aLl feOlO, and references therein). Most rel- 
evant to ou r work are the recent results pointing to a massive 
mol ecular dFischer et al.ll20 l ot IFeru glio et al. 2010) and neu- 
tral dRupke& Veilleuxll2011l) outflow in this galaxy. It is ar- 
gued strongly in these papers that radiative driving, rather than 
star formation or jet activity, power the observed outflows. In 
this paper we try to make a similar case for radiation-pressure 
driving, except that we focus on a radio-quiet obscured quasar, 
for which there is little evidence of massive star formation in 
the optical spectrum. By focusing on obscured, radio-quiet 
systems, we hope to isolate unambiguously the sources of 
acceleration in the observed outflows. We also hope to test 
the hypothesis that AGN feedback occurs most vigorously in 
the o bscured phase (e.g.. ISanders et al.ll!988h iHopkins et al.1 
12003) . 

Since 2007, we have been conducting a study of the mor- 
phologies and kine matics of narrow-line regions in obscure d 
radio-quiet quasars (Zakamska et al. 20031 iReyes et al.! 2008). 
in search of tell-tale signs of the interactions between the 
supermassive BH in its radiative (rather than jet-launching) 
phase and its large-scale gaseous environment. In our pro- 
gram conducted with the Ma gellan telescopes dGreene et al.l 
120091 hereafter Paper I, and lGreene et al.ll20lU hereafter Pa- 
per II), we obtained 15 long-slit spectra of obscured quasars. 
We find that in the quasar luminosity regime probed by our 
observations (L[omj~ 10 42 - 10 43 erg s" 1 or inferred Lt>oi~ 
10 45 — 10 46 erg s" 1 ), the quasar easily ionizes the galaxy 
along any direction not obscured by dust. Furthermore, 
the gas is disturbed on galaxy-wide scales, with both broad 
emission-line widths and kinematics clearly inconsistent with 
galaxy rotation. Among these 15 sources, one object, SDSS 
J135646.10+102609.0 (hereafter SDSS J1356+1026), shows 
very extended (> 50 kpc) ionized line emission which is 
highly kinematically organized on ~ 20 kpc scales. In this pa- 
per, we focus on this source, and build a case for AGN-driven 
feedback using an analysis of the kinematics and geometry of 
the ionized gas. 



2. OBSERVATIONS AND BASIC OBJECT PROPERTIES 

Despite its anonymous telephone number, the enigmatic ob- 
ject SDSS J1356+1026 has already garnered quite a bit of in- 
terest in the recent literature. It was targeted for sp ectroscopy 
by the Sloa n Digital Sky Survey dYork et al.ll20 00) in the Data 
Release 4 ( Adelman-McCarth v et alj|2006l) . Due to its large 
[O III]A5007A (hereafter [O III]) luminosity and line-ratios 
indicative of an active nucleus, it was included in the cat- 
alog o f obscured (type 2) quasars published by IReyes etall 
(2008). The kinematic structure of the [O III] line (as seen 
within the 3" SDSS spectroscopic fiber) is complex , show- 
ing at least two components. As a result, iLiu et al.l (120101) 
included SDSS J1356+1026 in a sample of type 2 AGNs with 
two double-peaked emission lines. 
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FIG. 1. — The broad-band spectral energy distribution of 
SDSS J1356+1026. Error bars represent the data presented in Table 
1 , while the blue line simply connects the points to guide the eye. 

Such signatures can be due either to the kinematic struc- 
ture of the narrow-line region itself or to two quasars, each 
with its own narrow-line region, with a velocity offset. The 
galaxy morphology is clearly disturbed, suggesting an on- 
going merger (Figure |2] Ze/f). The presence of two merging 
galacti c nuclei in SPSS J1356+1026 is apparen t both in opti- 
cal dGreene et al.l 120091 1201 It iFu et al.l 1201 lal) and in near- 
infrared (IShen et al.l 1201 ll) observations, and strong high- 
ionization emission lines are associated with both (Paper II 
and this paper). 

The spectroscopic observations of SDSS J1356+1026 were 
obtained on UT March 18, 20 07 using the Low-Dispersion 
Survey Spectrograph (LDSS3; Alling ton-Smith et al.lfl994h 
with a 1" x 4' slit at the Magellan/Clay telescope on Las Cam- 
panas. The seeing was ~ 1". We observed with two slit posi- 
tions for one half hour at each position. Slit position 1 (PA1 
hereafter) was oriented north-south, while slit position 2 (PA2 
hereafter) was oriented at 45° to the west of PA1 (Figure[2} ■ 
We used the VPH-Blue grism in the reddest setting, resulting 
in wavelength coverage of 4300 - 7050 A and velocity reso- 
lution of (Jinst w 67 km s" 1 . In addition to the primary science 
targets, at least two flux calibrator stars and a number of ve- 
locity template stars were also observed. 

Cosmic-ray removal was performed usin g the spectroscopic 
version of LACosmic (Ivan Dok kum 2001). Bias subtraction, 
flat-field correction, wavelength calibration, pattern-noise re- 
moval (see Paper I), and rectification were performed using 
the Carnegie Observatories reduction package COSMOS3 
The flux calibration correction is determined from the ex- 
tracted stan dard star using IDL rou tines following methods 
described in Mathes on et al.l {2008) and then applied in two 
dimensions. In Paper I we demonstrate that our flux calibra- 
tion is reliable at the ~ 40% level. 

Polarimetric observations were made at the 1 .54m Kuiper 
telescope with SPOL spectropolarimeter (Schmi dt et al.l 

1 http://obs.carnegiescience.edu/Code/cosmos 
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FIG. 2. — Left: r— band Magellan acquisition image of SDSS J1356+1026 and its immediate environment. Red scale bar indicates 10" (22 kpc) and blue lines 
indicate positions of the two spectroscopic slits. N is up and E is to the left. The observation taken with the N-S slit is PA1 and the one with the SE-NW slit is 
PA2. Right: The two-dimensional spectrum centered on [O m] A5007 taken at the N-S slit position. The vertical scale is the same as that of the image on the left. 
The full horizontal span is 3842 km s -1 centered at 5010A (rest frame). The "bubble" feature that is the focus of this paper is the round feature to the south of the 
continuum sources. There is also a tidal feature (TT) connecting the merging galaxies with the companion to the north, with the extent demarcated by the vertical 
white line. 



1992b) on March 29, 2011, with the 3"-wide slit oriented in 
the E-W direction and an integration time of 32 minutes. The 
same instrument was again used at the 2.3m Bok telescope on 
May 27, 2011 with the 3"-wide slit oriented in the N-S di- 
rection and an integration time of 40 minutes. In both cases, 
the polarization measurements span 4000-7550A with a res- 
olution of ~ 20 A. The degr ee of linear polarization ( P) is 
corrected for statistical bias (War dle & Kronberg||1974l) and 
the position angle of the polarization (9) was calibrated by 
observing the interste llar polarization stand ard stars Hiltner 
960 and VI Cyg #12 dSchmidt et al.lll992al) . Further details 
of the general observation and red uction procedures for the 
spectropolarimetry can be found in Sm ith et al.l (|2003). 

We have compiled the spectral energy distribution of 
SDSS J1356+1026 based on publicly available catalogs, from 
radio to ultra-violet wavelengths (see Table 1; Figure [TJ. X- 
ray observations with both XMM-Newton and Chandra are 
currently underway. 

3. ANALYSIS OF OBSERVATIONS 

3.1. Observed emission 

We first introduce the basic observational components of 
SDSS J1356+1026 and its immediate environment, as seen in 
Figure |2] We observe two nuclei, the northern (N) nucleus 
and the southern (S) nucleus, separated by ~ 2.5 kpc (1.1") 
on the sky. We find distinct stellar continuum and ionized gas 
emission from each nucleus in SDSS J1356+1026. Although 
the continua are spatially resolved, the southern source is faint 
enough that we are unable to determine the velocity separation 
between the nuclei from the stellar spectra. To the north is a 
companion galaxy located 62 kpc (in projection on the plane 
of the sky) from the N nucleus. In velocity space, the stellar 
features of the companion galaxy are blueshifted by ~ 250 
km s _I with respect to the stellar features of the N nucleus. 



In addition to stellar continua, there is strong, spatially re- 
solved [O III] AA4959,5007 line emission, as well as H/3 and 
other high-ionization forbidden lines. We focus here predom- 
inantly on the [O III]. A ~ 51 kpc-long [O III] emission line 
feature (denoted 'TT' for 'tidal tail'; Fig. |2]l extends between 
SDSS J1356+1026 and the northern companion galaxy and 
somewhat beyond it. No associated stellar continuum is de- 
tected for this feature, nor is it seen in the broad-band image. 
Additional likely tidal features are seen in the broad-band im- 
age, such as the low surface brightness feature ~ 10" (22 kpc) 
SE of the nuclei. This feature is tentatively detected in [O III] 
emission in the two-dimensional spectrum (see feature Bv2 in 
Figure |3J. 

The focus of this paper is on the well-organized kine- 
matic feature to the South of SDSS J1356+1026, "the bub- 
ble" (see Figure [3). The feature appears as a ring in the two- 
dimensional spectrum, with a spatial extent of ~ 12 kpc and 
apparent velocity extent of ~ 460 km s -1 . The central veloc- 
ity of the bubble at its base is close to that of the N nucleus, 
~ 8 kpc (projected) from its center. Equidistant to the north 
of the N nucleus, we find three clumps spanning a comparable 
velocity range. The bubble and the clumps are pure emission- 
line features, with no stellar continuum detected in the spec- 
tra. Emerging at a 45° angle to the bubble are two 'wisps' of 
emission associated with the same outflow. 

We subtract the stellar continuum and co-add the spectrum 
in velocity space to evaluate the [O III]/H/3 line ratio as a func- 
tion of position along the slit (Figure|4]i. The ratio is constant 
along the slit, with [O III]/H/3 ~ 10, suggesting that the gas is 
photoionized by an accreting BH. This line ratio is constant 
for a remarkable number of components, including both the 
S and N nuclei, bubble, and the northern clumps. We have 
not corrected these observations for extinction , but note that 
assuming Ry = 3.1 and the reddening law of ICardelli et alj 
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TABLE 1 
Spectral Energy Distribution 
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NOTE. — Spectral energy distribution of SDSS J1356+1026 from archival data. SDSS fluxes are the average of Petrosian and cmodel fluxes. SDSS errors are half of 
the difference between Petrosian and cmodel fluxes combined in quadrature with the nominal error on the Petrosian fluxes. The large uncertainty in the i band is presumably 
due to the highly disturbed morphology of the object, so that the standard SDSS model fits fail. FIRST and NVSS fluxes are 59.58 mJy and 62.9 mJy, correspondingly. 
The table entry lists the average of the two. with the uncertainty being half the difference between the two. Although the nominal threshold for the 74 MHz catalog is 
5cr = 0.70 Jy, there is no hint of any flux at the source position in the survey stamps, so a 2a upper limit is adopted. The far-infrared luminosity of SDSS Jl 356+1 026 
as estimated using the fitting formula by Sanders & Mirabel 1 1996) is 2.68 X 10 45 erg/sec (the uncertainty is dominated by IRAS flux calibration, ~ 20%). We use an 
h = 0.7, n,„ = 0.3. SI a = 0.7 cosmology. 



( 1989), there is only an 8% difference in the extinction of H/3 
and [O III]. Thus, if we invoke (for instance) higher levels of 
extinction in the galaxy center, we would still only expect the 
[O III]/H/3 ratio to decrease by < 10% in the central region. 

As argued in Paper II, the most plausible explanation for the 
constant line ratios is that both nuclei host obscured quasars 
with similar luminosities. Another possible scenario is that 
only one quasar is present and it photo-ionizes the other nu- 
cleus; however, we find that the fine-tuning required to sat- 
isfy the available observables makes this hypothesis unlikely 
(Paper II). Future observations (in X-rays with Chandra or in 
the radio with ~1" resolution) will help distinguish between 
these two p ossibilitie s. 

Recently iFu et alj d201 lbl) presented integral field obser- 
vations of the 6"x6" field around SDSS J1356+1026 (e.g., 
largely excluding the bubble itself). It is not obvious from 
their [O III] map that they see emission line sources associ- 
ated with both nucleus N and nucleus S. They conclude that 
SDSS J1356+1026 harbors only a single active black hole in 
nucleus N, and that the multiple kinematic components are 
due to large-scale motions of the gas. It is possible that our 
slit did not intercept the peak emission from nucleus N, thus 
explaining why the continuum strengths of nucleus S and nu- 
cleus N are roughly comparable in our spectrum, even though 
source N is considerably more luminous (Paper I). We still 
believe we are seeing a distinct [O III] peak associated with 
nucleus S, but cannot perform more detailed comparisons at 
this time. 

Taking the spectrum of the bubble alone and subtracting 
the continuum, we use the velocity profile of the [O III] 
line to fit all other emission features in the spectrum and 
thus calculate line ratios without modeling the complex ve- 
locity structure of the lines. In particular, we find a ra- 
tio of [O III] AA4959,5007/[O III] A4363= 90.5 ± 7.2 and 
therefore derive an electron temperature T e = 13,500 ±500 



K (lOsterbrock & Fer land 2006. Furthermore, we detect the 
[Ar IV] AA471 1,4730 doublet with a rati o of 1.32 ±0.37, con - 
sistent with its low-density limit of 1.4 dKrueger et al.ll 1 970b . 
Because the critical density for this transition is fairly high, 
even if we take our measurement as a lower limit on the ratio, 
it places only a weak limit on the electron density (n e < 850 
cm -3 ). The ratios of lines of the Balmer series are consistent 
with no reddening. 

The Magellan spectrum does not cover the density diag- 
nostics [Oil] A A 3726,3729 and [S II] A A 6716,6731, but both 
of these features are covered by the SDSS spectrum taken in 
the 3"circular aperture centered on the nuclei. The compo- 
nents of the [O II] doublet are separated only by 240 km s~' 
and are far too blended to determine their ratio, given that 
FWHM [0 iii]= 652 km s" 1 in this aperture. The [S II] doublet 
has an acceptable velocity separation, but the velocity struc- 
ture of the [O III] emission line does not provide a good fit to 
this feature, so we cannot determine the line ratio in the same 
model-independent way as for the Magellan spectrum. We 
measure the peak ratio of the two components of the doublet 
to be 1.20, which is likely to be a lower limit on the actual ra- 
tio because of blending. Taking [S II] 6716/6731 >1.2 yields 
n e < 250 cm -3 . 

The object is strongly linearly polarized, with the polariza- 
tion fraction reaching 8% at shortest observed wavelengths 
(Figure |5}. These observations are consistent with the hy- 
pothesis that SDSS J1356+1026 is powered by a luminous 
quasar whose emission is largely blocked along the line of 
sight to the observer, but escapes along some other directions, 
and the n scatters off the i nterstellar material wi thin the host 
galaxy (ISmifh et al.ll2003t IZakamska et a l. 2005) and reaches 
the observer. In particular, the emission from the quasar es- 
capes predominantly in the N-S direction in the plane of the 
sky, as indicated by the ^90° PA of the plane of polariza- 
tion. Furthermore, the wavelength dependence of the polar- 



SDSSJ1356+1026 



5 




5 -• 



,,,,,,, 

- ^ Bubble x 2 ^J\^ 


{ " 


S x 0.5 y \ _/ 


-X - 






Clumps x 50 nn 


Companion x 50 fl n 

_, , 1 , , , , 1 , , i , 1 , , , , 1 





4850 4900 4950 5000 5050 
Rest Wavelength (A) 









i "3rl f i |T 










s 


ill 







5 - 



2.5 




4850 4900 4950 5000 5050 
Rest Wavelength (A) 



FIG. 3. — Two-dimensional spectra (left) and one-dimensional extractions along the labeled regions (right) for slit position PA1 (top) and PA2 (bottom). The 
two-dimensional spectra are centered on 5007 A and span 80A (rest-frame) in the spectral direction. Each image is 23 " (50 kpc) in height, although the spectrum 
along PA2 extends 15 kpc further to the south to highlight Bvl and Bv2 (the 'wisps'). In the two-dimensional spectra we label features of interest, including, from 
PA1 (top left) - B: bubble; S: southern nucleus; N: northern nucleus; C: clumps. We also show the one-dimensional spectrum of the northern Companion seen 
in Fig. |5]with a 65 kpc projected distance from our program object. We label in PA2 (bottom left) - Bvl: wisp feature 1, Bv2: wisp feature 2, and C: the small 
amount of flux in the same region as the clumps. The one-dimensional spectra (right) are extracted in the windows shown in the two-dimensional image, and the 
labels correspond, except that the B spectrum bottom right is the sum of Bvl and Bv2. The spectra have been multiplied by the factors indicated on the plot for 
ease of display. In addition, for the PA2 spectra (bottom right), we overplot the corresponding features from PA1 to demonstrate the relative velocities between 
the components. For instance, note that the bubble and the wisps Bvl and Bv2 have a similar range of velocities. The artifact in the Clump and Companion 
spectra along PA1 are due to a gap between the two LDSS3 CCDs. 



ization fraction can be attributed to the dilution of the scat- 
tered quasar light by the unpolarized emission from the ex- 
tended narrow-line region and starlight from the host galaxy. 
Indeed, the dilution of the polarized flux by starlight increases 
to the red, as expected from a red stellar population, and an 
abrupt minimum in P is seen around [O III] (Figure |5). The 
generally higher values of P for the N-S slit orientation of the 
May 27, 201 1 observation are likely caused by the fact that 
the spectral extraction aperture runs along the axis of the pre- 
sumed scattering cone, thereby including a larger fraction of 
the scattering regions than the earlier measurement made with 
the Kuiper Telescope. Other small discrepancies between the 
two spectropolarimetric observations are likely caused by dif- 
fering positions of the object within the slit and by the differ- 
ence in the slit direction, suggesting that polarized structures 
within the source are partially resolved. 

3.2. Observed Kinematics and Geometry 



We propose that SDSS J1356+1026 is powering a large- 
scale outflow, which is entraining the warm ionized gas that 
we see radiating in [O III]. Here we present the observed kine- 
matics along each slit position that we will use to estimate the 
energetics of the outflow. 

Starting with the bubble along the first (NS) slit position 
PA1, we take two-pixel-wide extractions along the spatial di- 
rection starting from nucleus S and working our way to the 
bottom of the bubble (all within the bubble as demarcated in 
Figure|3]l. We use Gaussian fits to identify the primary veloc- 
ity components at each spatial position. These velocities are 
plotted in Figure [6] ignoring some weak components. This 
velocity pattern, increasing symmetrically to a projected ve- 
locity of ~ 250 km s -1 , is the signature of an expanding sphere 
of gas with limb brightening. 

In addition to the bubble, there are three clumps to the north 
of nucleus N that we see in slit position PAL We note that (a) 
their central velocity is similar to that of the bubble and (b) the 
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maximum velocity extent across the clumps is comparable to 
that seen in the bubble (see Figure |3). Therefore, we suggest 
that these clumps are the northern side of the same outflow. 
The bubble is systematically blueshifted by ~ 215 km s" 1 rel- 
ative to the clumps; therefore, the axis of symmetry of the 
outflow is likely not quite in the plane of the sky but is tipped 
so that the bubble is expanding slightly towards us and the-£ 
clumps are from the component moving slightly away from p 
us (see Figure |5). The northern and southern components of g 
the outflow are discontinuous; we suggest three possible cor- 
responding geometries. One possibility is that we are seeing 
two unrelated components, but given the similarity of the sys- 
temic velocities and the velocity spread we find this unlikely. 
Another possibility is that our slit did not pass through the 
center of the outflow, but instead was 1-2" (2.3-4.6 kpc) off- 
set (see arguments above suggesting that the PA1 slit was not 
centered on nucleus N). In this case an outflow that forms a 
figure '8', with the center taken out, could appear as a pair of 
disconnected features. Yet another possibility (that we favor 
on the basis of kinematic model presented in the next section) 
is that the merger provides plentiful gas in the inner several 
kps of the system, pinching the flow, which finally breaks free 
of this dense gas on either side, as shown in Figure 5. 

The signal-to-noise ratio along PA2 is not very high, and 
thus we do not measure the velocities as a function of position 
in this case. The two linear structures to the south of nucleus 
N, Bvl and Bv2 (Figure show a clear trend of increasing 
line-of-sight velocity with projected distance. These wisps 
represent gas that is escaping at a 45° angle to the gas in the 
bubble. 

3.3. Kinematic Model 

We see outflowing gas to the north, south, and south-east 
of the nuclei of SDSS J1356+1026. The emission-line gas 
is clearly expanding away from the nuclei. Is it being ac- 
tively driven from the nuclei and accelerating away, or is it 
decelerating as it slams into the interstellar or intergalactic 
medium? The linearly expanding structures along PA2 (wisps 
Bvl and Bv2) are the tell-tale signs of an accelerating outflow 
(IVeilleux et al.lll994T) . Both a conical outflow with V3d oc r^o 
and a spherical outflow blown from one side with vt,d oc r\ D 
will produce v oc r when projected onto the line of sight (v) 
and on the plane of the sky (r). The cap of the bubble along 
the southern-most edge is produced by a shell where the out- 
flow interacts with gas outside the galaxy, either tidal debris 
or the intergalactic medium. Because of projection effects, the 
cap is brightest in the plane of the sky, where its line-of-sight 
velocity is zero, so we do not know whether the shell is still 
accelerating or whether it has already started to decelerate. 

If we assume that there is a single accelerating agent caus- 
ing all of the outflowing components, then a credible picture 
of the geometry can be built. Due to interactions with the am- 
bient medium of varying density, the outflow is complicated 
enough that we cannot uniquely determine its geometry from 
two slit positions, but we are able to bracket the outflow ve- 
locities using a reasonable range of models. Our primary goal 
is to explain the symmetric velocity structure along the bubble 
in a way that is also consistent with the clumps to the north 
and the linearly increasing velocity structures ('wisps') to the 
south-east along PA2. 

We first consider an accelerating (v3d oc rso) biconical out- 
flow, which is a natural geometry for an obscured quasar and 
easily explains the wisps. 
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FIG. 4. — The spatial profiles of [O in] and H/3 emission along the N-S slit. 
The [O m]/H/3 ratio is constant and close to 10. 

This model fails to reproduce the concave velocity structure 
(the southern "cap") observed in the bubble (Figure |6}. An- 
other straightforward model is a spherical outflow centered on 
nucleus N. If only the shell is emitting, it is then difficult to re- 
produce v oc r at small radii in the bubble and the wisps. Both 
the cap and the accelerating features are naturally reproduced 
in the framework of a bubble blown from its side, not from 
its center, perhaps due to the high density o f matter near the 
energy injection point ( Veille ux et al.l [19941) . This model is 
reminiscent of outflows in starburst galaxies, in which prop- 
agation in the plane of the galaxy is strongly suppressed and 
the outflows escape p redominantly perpendicular to the plane 
( Veill eux et al.ll2005l) . Note that if the outflow axis is inclined 
toward the plane of the sky, then the velocity offset between 
the clumps and the bubble (relative to the maximum span in 
the velocity of the bubble) leads to an estimate of the inclina- 
tion of the axis of ~ 20°, which could in principle boost the 
space velocities relative to those assumed here. 

A simple version of this model is shown in Figure [6] We 
assume that the flow is spherical, with the origin on the side 
of the sphere and the axis in the plane of the sky, that only 
the shell is emitting and that V30 oc r" D where n = 1,1.5,2. 
Then we project the system onto the sky and the velocities 
along the line of sight and plot the model curves along with 
the observations of the bubble in SDSS J1356+1026. We 
perform a \ 2 fit to the observed position-velocity diagram 
taking 1 < n < 2.5. We exclude all data points with veloc- 
ities greater than ±400 km s , since they are not treated 
by our model. We further use a minimum velocity error of 
20 km s" 1 (a fraction of a pixel). Models with 1.2 <n < 1.5 
all yield a reduced \ 2 differing from the minimum by less 
than one. The corresponding range of projected velocities is 
250 < V30 < 1000 km s . At the lowest velocity there is no 
projection effect: the maximum observed velocity is also the 
maximum space velocity. This is a strict lower limit to the 
energetics. 

We postulate that all of the line emission is part of a quasi- 
spherical outflow blown by the AGN, that extends from the 
southern-most edge of the bubble to the clumps in the north 
(illustrated schematically in Figure|6|. However, due to inter- 
ference from gas and dust in the central pair of galaxies, the 
outflow cannot escape freely in all directions. Along PA2, the 
wisps represent gas that has broken free and is in free expan- 
sion. In contrast, along the N-S direction (PA1) the outflow 
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runs into the two galaxies. The ring shape in the PA1 2D spec- 
trum results from gas breaking out of the high density mate- 
rial associated with galaxy S and accelerating away from the 
merger with V30 oc r]^. In addition, the structure is discon- 
nected from its origin (nucleus N) because of the intervening 
matter. The clumps in PA1 are the northern counterpart of the 
same flow, produced closer to nucleus N than the bubble since 
there is less material to penetrate. In this picture, the similar 
velocity spread seen along the bubble, clumps and wisps (Fig. 
are naturally explained as arising from the same outflow. 
The opening angle of the outflow is no less than the 45° angle 
between PA1 and PA2, since both slit orientations captured it. 

3.4. Energetics 

With our crude outflow model, and thus deprojected veloc- 
ities, we can now estimate the energy injection E required t o 
power the observed bubble, following Nesva dba et al.l (2006). 
One method is to estimate the kinetic energy in the shell, 
based on an estimate for the gas mass and the observed ve- 
locity, and then divide by the dynamica l time to estimate an 
energy injection rate (e.g.. lMartiri fl998). The other approach 
is to treat the bubble like a supernova remnant, and estimate 
the energy needed to expand the bubble into a low-density 
ambient medium. In both cases we are forced to make as- 
sumptions, but with the two methods we can at least bracket 
the kinetic luminosity of the outflow. 

First, we derive a lower limit on the energy in the outflow 
based on its apparent kinetic energy. The dominant uncer- 
tainty here is the mass of gas in the bubble. In principle, with 
an electron density, we can turn the observed H/3 luminosity 
(Lh/3 = 7.5 x 10 41 ergs s" 1 ) into a hydrogen gas mass. As ex- 
plained in detail in Paper II, we have two methods to estimate 
gas density in the extended nebulae of obscured quasars. One 
is the standard diagnostic line ratios, which give density esti- 
mates n e < 1000 cm" 3 . Unfortunately, this method is biased 
toward regions of the brightest line emission, and thus regions 
of the highest density. Furthermore, in clumpy media, density 
diagnostics often yield values close to the critical density in a 
wide range of conditions (Paper II). The other method is to use 
the observations of scattered light at large distances from the 
nucleus to infer a gas density, but this has only been possible 
for a couple of objects observed with the HST and with partic- 
ularly well-measured scattered light nebulae (Zakamska et al. 
2006). This method gives a mean electron density of n e < 1 
cm -3 . Clearly, the two estimates are discrepant, suggesting 
that gas clumping is significant. Defining the clumping factor 
as k = (n 2 ) I (n e ) 2 , we can write the mass in gas as: 



M g = 1.7 x 10 M Q -L^(n e y 



(1) 



In section l3TI we used diagnostic line ratios to find that in the 
nuclear region n e < 250 cm -3 and that [O III]/H/3 is constant 
throughout the nuclei and the bubble. Since at the location of 
the bubble the ionizing source is 5-10 kpc away, the density 
must drop with radius by at least a factor of 20-100 to main- 
tain the same ionization parameter (the ratio of the number of 
ionizing photons to electrons, U). This decrease can plausibly 
be somewhat offset by a clumping factor of similar magni- 
tude. Therefore, our final mass estimate is M g >5x 10 7 M Q . 
This mass is a lower limit, because it takes into account only 
the ionized gas that produces emission lines. Any lower den- 
sity gas at higher temperature filling the bubble would not be 
included in this estimate. 
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FIG. 5. — Top: Polarization fraction. Bottom: polarization position angle 
(degrees E of N) as a function of rest-frame wavelength. Black skeletal points 
are from the Kuiper telescope observations, with the 3"-wide slit in the E-W 
direction, and and grey open circles denote binned spectropolarimetry from 
the Bok telescope, with the 3"-wide slit oriented N-S. 

On the other hand, without a measurement of the clumping 
factor k, we cannot say for sure what fraction of the ob- 
served volume is filled by high-densit y gas. In the case of 
the observed bubble in NG C 3079 (e.g., Filippenk o & Sargentl 
[1991 lVeilleuxeFaT][T99ll) . the energetics estimates were re- 
vised downwards based o n high-resolution HST imaging of 
the ionized gas filaments (ICecil et a l. 2001). We plan to ob- 
tain integral-field observations of SDSS J1356+1026 in the 
Spring, which will hopefully help settle the outstanding am- 
biguities. 

Based on the extent of the bubble, R = 12 kpc, and the max- 
imum deprojected velocity of v = 250- 1000 km s" 1 , we find 
the dynamical time, fa ~ R/v » 10-40 Myr. It is interesting 
to note that this timescale is similar to ty pical AGN lifetimes 
(Ma rtini & Weinbergl20M lMartini 2004). With our previous 
estimate of the mass of the emitting gas, this corresponds to a 
mass outflow rate of M ~ M g /fa > 2-5 M Q yr" 1 . The corre- 
sponding kinetic energy is E^ n ~M g v 2 /2 > 3-50 x 10 55 ergs. 
For a constant energy injection rate over the lifetime of the 
bubble, we find E > 3- 10 x 10 41 ergs s" 1 . This is the kinetic 
energy of the emission line gas alone and, therefore, a strict 
lower limit on the total energy of the outflow. 

The energy injection rate required to expand the bubble into 
a low-density medium, no, in the galaxy halo can be calculated 
in a manner simi lar to the Sedoy-Taylo r solution for super- 
novae. Following iNesvadba et al. (2006), we assume that the 
radiative losses are minimal and that energy injection rate is 
constant rather than a single explosion. We then find 



E « 1.5 x (n/4TT)l0 46 R 2 w v 3 mo n Q ergs s" 1 



(2) 



where the extent of the bubble Rio is in units of 10 kpc, viooo is 
the expansion velocity in units of 1000 km s" 1 and SI /4tt is the 
covering factor of the outflow in steradians. Taking Rw=l .2, 
no = 0.5 cm" 3 , and viooo = hE » (0/4tt) x 1.1 x 10 46 erg s _1 . 
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If we instead adopt the lower limit on velocity (i.e., no depro- 
jection)of vjooo = 0.25, we find£ w (Q/4tt) x 2 x lO^erg s" 1 . 

In addition to the deprojected velocity, the major uncer- 
tainty in this value comes from the lack of good measurements 
of «o and f2. Both slit positio ns, separated b y 45°, captured 
the signature of the flow and Fu et al. (201 fa|) find a strong 
component of [O III] emission to the northeast of N, which is 
not covered by our observations. We therefore conclude that 
£7 is not too small. It is much harde r to estimate the proper 
value for «o- The HST observations of Zakamska et al. ( 2008) 
alluded to above find clear evidence for electron scattering off 
of a medium with n e < 1 cm" 3 , even at these large projected 
distances of ~ 10 kpc from the nucleus. We clearly see evi- 
dence for denser gas from the observed line emission. Thus, 
we view our choice of «o to be a conservative lower limit. 

The kinetic energy injection rate of the outflow is bracketed 
by the values 10 42 erg s" 1 (where we likely underestimated 
the amount of mass by a large factor) and 10 46 erg s" 1 (where 
we likely overestimated the outflow covering factor and the 
density of the ambient medium by smaller factors). We will 
therefore adopt a fiducial range of 10 44 " 45 erg s" 1 for the ki- 
netic energy injection rate in the discussion that follows. It is 
interesting to note that our upper limit on the kinetic energy 
injection rate is roughly consistent with our estimate for the 
radiative bolometric luminosity of Lboi ~ 10 46 ergs s . 

3.5. Can the outflow be driven by star formation? 

Can the energy of the outflow be plausibly supplied 
by any source other than accretion onto the supermassive 
black hole? Vigorous star forma tion can lead to large- 
scale outflows ( Veil leux et al.ll2005[) that are observable both 
in ionized gas in em ission (e.g.. IChevalier & deed fT98l 
Heckman et al.lll990j) and in neut ral gas in absorption (e.g., 
Rupkeeta In an active merger like 

SDSS J1356+1026, with plenty of gas available (as evidenced 
by our emission line observations), it is natural to expect 
that some on-going star formation is present. While from 
the high-ionization line ratios (such as the uniformly high 
[O III]/H/3 ~ 10) it is clear that the photo-ionization of the gas 
is completely dominated by the quasar, star formation may 
still be the primary driver of the kinetic outflow. Unfortu- 
nately, despite the availability of abundant multi-wavelength 
data and deep spectra of the host galaxy, the relative contri- 
butions of star formation and obscured quasar activity to the 
bolometric output of SDSS J1356+1026 are very hard to de- 
termine. This is a generic difficulty arising due to the fact that 
most of the flux from the object is coming out at mid- and 
far-infrared wavelengths, where the emission is due to ther- 
mal radiation of dust particles heated by young stars or the 
quasar in a manner notoriously insensitive to the spectrum of 
the underlying power source. 

The bolometric luminosity of the quasar, as expected 
from its [O III] luminosity dLiu et all 120091) . is 1.5 x 10 46 
erg s" 1 , although such scalings have dispersion of > 0.5 
dex. The bolometric luminosity estimated using the 25 
and 60/xm IRAS fluxes and quasar bolometric corrections 
dRichards et al.ll2006l) ranges between 7.2 x 10 45 -2.4 x 10 46 
erg s" 1 . The reasonable agreement between these values indi- 
cates that most, if not all, of the mid- and far-infrared lumi- 
nosity from SDSS J1356+1026 (2.7 x 10 45 erg s" 1 ) can plau- 
sibly be produced by the quasar activity rather than star for- 
mation. The high degree of polarization likewise suggests that 
the AGN is dominant. 



If, however, all of the infrared luminosity of this source 
were powered by star formation, it could just barely power 
the observed outflow. We use the calculation of the me- 
chanical ene rgy injection rate from a co nsta nt rate of star 
forma tion bv lLeitherer & Heckmanl (11995b and lVeilleux et al.l 
(2005). For every solar mass of stars formed, the combi- 
nation of stellar winds and supernovae can produce at most 
^mech ~ 7 x 10 41 ergs s" 1 at time t rs 10 7 yr. Thus, to produce 
a 10 44 " 45 erg s" 1 outflow would require > 100- 1000M Q yr" 1 
of star formation and predict a star-formation luminosity of 
i-FiR = 3 x 10 45 " 46 ergs s . Thus, star formation is marginally 
consistent with powering the outflow at the lower limit of the 
energy range. On the other hand, a starburst of > 100 M Q yr" 1 
is hard to reconcile with our optical spectroscopic and spec- 
tropolarimetric observations. We see no evidence in the spec- 
trum for an intermediate age population (i.e., the Balmer ab- 
sorption lines are consistent with an old stellar population), 
and the blue contin uum that we d o see is consistent with 
scattered light (e.g.. lLiu et alj |2009). Therefore, we consider 
power from star formation as unlikely given the currently 
available data. 

3.6. Can the outflow be jet-driven? 

Radio flux densities at 1 .4 GHz can be obtained from the 
Faint Images of th e Radio Sky at Twenty cm survey (FIRST; 
[ Becker et alj|"l9 95) and the NRAO VLA Sky Survey (NVSS; 
Cond onet al.ll 1998b . The FIRST detection is unresolved - the 
nominal deconvolved siz e of 1.2x 0.6" i s smaller than half 
of the beam size of 5.4" ( iHo & UlvestadlHOOl) . The FIRST 
peak flux is 57.90 mJy beam" 1 and the integrated flux is 59.58 
mJy; the ratio of these two measures is consistent with that 
seen in point sources (Ivezi c et alj 120021; iKimball & Ivezid 
2008). The high-resolution FIRST observations tend to re- 
solve out the low surface-brightness extended emission, so the 
comparison with lower-resolution NVSS observations (45" 
beam) provides yet another measure of extended emission. 
The NVSS catalog flux is 62.9 mJy, and again the difference 
between NVSS and FIRST- measured fluxes is consis tent with 
that seen in point sources (Kimball & I veziq 120081) . From 
these observations, it appears that SDSS J1356+1026 lacks 
extended radio emission on > 4 kpc (2") scales at the few 
mJy level. 

U sing the radio-infra red correlation of star-forming galax- 
ies dHelou et alj I1985[) a n d its convenient normalization 
presented by Mori c et all ([2010), we find that no more 
than 20% of the observed radio emission can be due to 
star formation. Thus, most of the radio emission is at- 
tributable to one or both quasar nuclei. Nevertheless, 
the positi on of the SPSS J1356+1026 in the P O III] -radio 
diagr am dXuetalJ [19991: iZakamska et all l2004t lLal & Hoi 
2010) puts it squarely into the ra dio-quiet regime , 1.5-3 
dex below radio-loud sources in iNesvadba et al] (2008) 
and iFu & Stockton] d2009l) . For direct comparison with 
lXuetal.1 (11999). we use a radio spectral index of a = 
-0.47 (defined as F v oc v a , see below) to convert a rest- 
frame-corrected value of log(fi.4GHz, erg s" 1 Hz" 1 ) = 31.32 
into log(P 5G Hz, e rgs" 1 Hz -1 ) = 31.0 6 at log(L [0 ni], erg s" 1 ) = 
42.78 taken from Re yes et ail (120081) . Obscured quasars with 
such radio luminositi es have only a ~ 10% chance to have ra- 
dio extents > 3 kpc dLal & Holl2010l) . consistent with the lack 
of any sign of extended emission in SDSS J1356+1026. In 
principle, we cannot rule out a < 4 kpc-scale jet as the acceler- 
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FIG. 6. — We show our simple kinematic model (left) and an accompanying cartoon of the outflow region (right). Left: Position-velocity diagram centered 
on the the [O III] bubble, as demarcated in Figure [3] Filled black circles indicate velocity peaks measured at intervals of 0'.'4. We assume that the accelerating 
agent is offset from nucleus N by ~ 2", thus the outflow is detached from the continuum source. While the acceleration arises from nucleus N, we show only the 
outflowing region, starting at the open circle at the origin. Overplotted is a schematic model of a quasi-spherical outflow that is decelerated at its base as it passed 
through continuum sources N and S. The velocity along the flow increases with radius from source N with a power-law dependence v oc r". We plot models 
with n = 1 (red dash-dot), n = 1 .5 (blue dotted) and n = 2 (green dashed). This model yields a maximum outflow velocity of ~ 1000 km s" 1 . Right: A schematic 
cartoon illustrating our best guess of the system geometry, including the two continuum sources, the Bubble to the south, Clumps to the north, and the slight tilt 
in the plane of the sky relative to the observer. The circle indicated the origin of our model to the left. 

ionized by the quasar along an unobscured direction. High 
optical linear polarization is observed, also indicating that the 
emission-line regions are powered by the AGN. In addition, 
we observe a pair of expanding shells, or bubbles, centered 
on the northern nucleus and expanding into the intergalac- 
tic medium. The remarkable scale of the expanding structure 
(~ 40 kpc in total extent), its high expansion velocities (ob- 
served velocity range of 460 km s -1 , with physical velocities 
likely reaching 1000 km s" 1 ), and its large opening angle and 
luminosity imply that the kinetic energy of the outflow likely 
exceeds > 10 44 " 45 erg s4 

We find that while star formation may be occurring concur- 
rently with quasar activity in SDSS J1356+1026, it is energet- 
ically insufficient to power the observed outflow, likely by an 
order of magnitude, as further supported by the polarization 
observations. The merging quasars are radio quiet, and there 
is no evidence for extended radio emission; therefore, the out- 
flow is unlikely to be powered by a relativistic jet in a manner 
similar to outflows seen in some radio galaxies. We propose 
that SDSS J1356+1026 may be an example of on-going feed- 
back from a quasar in its most common radiative phase. Such 
feedback has been postulated to acc ount for the colors an d lu- 
minosities of massive galaxies (e.g., Springel et al. 20051), but 
it has been difficult to assemble direct evidence of interactions 
between radio-quiet quasars and their large-scale gaseous en- 
vironments aside from individual cases noted in the introduc- 
tion. SDSS J1356+1026 provides compelling evidence that 
accretion energy can drive wide-angle outflows of gas beyond 
the confines of the host galaxy. 
The presence of extended ionized emission by itself is not a 



ating agent, as has been proposed for NGC 3079 dCecil et alj 
2001). On the other hand, the outflow we observe is consid- 
erably more extended and more energetic than that in NGC 
3079, so it is not clear that energy injection on 1 kpc scales 
can power the observed 20 kpc outflow. Deeper radio inter- 
ferometry at 1" resolution would settle the issue. 

The source is und etected in the Texas survey at 365 MHz 
(Doug las et al.lll996|) and in th e VLA Low-Frequency Survey 
at 74 MHz dCohen et ai]|2007l) . with the 74 MHz observation 
providing a stronger upper limit, so we focus on this obser- 
vation. This implies a limit on the spectral index a > -0.47 
(hence the use of this value in the calculation above). Such 
a flat spectrum effectively ex cludes a "relic" radio source 
dKom issarov & Gubanov 1994])- 

High-resolution, high-sensitivity radio observations would 
be useful in establishing the presence of two quasars and of 
any extended radio emission, especially on 5-10" scales. 
Since the separation of the nuclei is ~ 1", observations with 
VLA would be well-suited for these purposes. The archival 
radio observations available to us at the moment indicate that 
the outflow in SDSS J1356+1026 is unlikely to be jet-driven. 

4. CONCLUSIONS 

SDSS J1356+1026 is an on-going galaxy merger at red- 
shift z = 0.123 that contains at least one obscured quasar in 
its northern nucleus. The bolometric luminosity of the object, 
estimated at 1.5 x 10 46 erg s " 1 , places it among th e more lu- 
minous nearby quasars (e.g., Hop kins et al.| [2007). We have 
conducted a study of the kinematics and spatial distribution 
of the ionized gas in this source. We find a ~ 50 kpc tidal fea- 
ture close to the systemic velocity of the object that is photo- 



10 



GREENE, ET AL. 



signature of quasar feedback (Sto ckton & MacKenty||1987h . 
If one could obtain arbitrarily sensitive observations, one 
would be bound to find small amounts of gas, and if this 
gas finds itself along an unobscured direction to the quasar 
it can be photo-ionized o ut to large distances. For example, 
iVillar-MartinetaLl d2010) find a 180 kpc emission-line fea- 
ture near an SDSS type 2 quasar, but its spatial distribution 
and kinematics make it clear that the feature is tidal debris il- 
luminated by the quasar. Images of this source t aken with the 
Hubble Space Telescope (Zakamska et al. 2006) demonstrate 
that the scattering regions (and therefore illumination cones) 
have the same E-W orientation as the 180 kpc tidal tail. We 
find a very similar feature in SDSS J1356+1026 extending in 
the N-S direction, which also happens to be the direction of 
quasar illumination, as indicated by the polarimetric observa- 
tions. Rather than swept from the galaxy by quasar feedback 
to 50-100 kpc away, the gas seen in these two cases is most 
likely tidal debris which happened to be along an unobscured 

line of sight to the qua sar. 

Similarly, [Liu et al.l (120091) find several obscured quasars 
with companions out to > 30 kpc away from the quasar show- 
ing high equivalent width emission lines with high-ionization 
line ratios. Again, rather than instances of quasar feedback, 
these occurrences are probably small companion galaxies 
which happen to fall within one of the quasar illumination 
cones. 

What sort of signatures specifically implicate quasar feed- 
back? The absence of stellar continuum is a necessary (al- 
though not a sufficient) requirement, as it helps rule out illu- 
minated companion galaxies and tidal features. Gas velocities 
that exceed the escape velocity of the galaxy are a clear sign, 
but since the outflows from type 2 quasars are directed in the 
plane of the sky, the chances of observing such high veloc- 



ities in these objects are low. The presence of a kinemati- 
cally organized and symmetric structure, such as that seen in 
SDSS J1356+1026, is a very compelling clue, as it completely 
rules out an instance of a gas inflow illuminated by the quasar 
(while the nature of radial velocity observations technically 
makes inflow and outflow scenarios degenerate, it is hard to 
imagine the gas conspiring to accrete in such symmetric man- 
ner). The high linear-polarization fraction of P = 8% pro- 
vides unambiguous evidence that there is a luminous obscured 
quasar in SDSS J1356+1026 and severely limits the fraction 
of the bolometric luminosity of this source that may be due 
to star formation. Furthermore, through the measurement of 
the polarization position angle, polarimetric observations in- 
dicate that the outflow direction (N-S) is well aligned with 
the predominant direction of quasar illumination. Clearly a 
narrow-band image or integral-field spectroscopy would pro- 
vide a more complete picture of the energetics. 

We did not find similarly well-organized kin ematics in the 
rest of the Magellan sample dGreene et al.l201 il) . but the large 
velocity widths and the lack of rotation seen even in promi- 
nent disk galaxies observed along their semi-major axes make 
quasar feedback the likely explanation in those cases as well. 
Our work to study three-dimensional outflows in luminous 
obscured quasars is ongoing with integral-field spectroscopy, 
which we hope will provide better statistics on the population 
as a whole (Zakamska, N. L. et al. in preparation). 

The referee made valuable comments that significantly im- 
proved this manuscript. We gratefully acknowledge useful 
conversations with Daniel Proga. We also thank L. C. Ho and 
A. J. Barth for invaluable assistance in the earlier stages of 
this project. RS.S. acknowledges support from NASA/Fermi 
Guest Investigator Program grant NNX09AU10G. 
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